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[215] Attorney Docket No.: UNI-107 Patent 

5 ASYMMETRIC RADIO FREQUENCY COILS FOR MAGNETIC RESONANCE 

FIELD OF THE INVENTION 

This invention relates to radio frequency coils for magnetic resonance (MR) 
applications. In particular, the invention is directed to asymmetric radio frequency 
10 coils for magnetic resonance imaging (MRI) machines. 

In certain of its aspects, the invention provides methods for designing radio 
frequency coils for magnetic resonance applications which may be symmetric or 
asymmetric. 

The radio frequency coils of the invention may be used for transmitting a radio 
1 5 frequency field, receiving a magnetic resonance signal, or both transmitting a radio 

frequency field and receiving a magnetic resonance signal. When the radio frequency 
coil serves a transmitting function, it will normally be combined with a shield to reduce 
magnetic interference with external components of the magnetic resonance imaging 
system. 

20 BACKGROUND OF THE INVENTION 

In magnetic resonance imaging (MRI) applications, a patient is placed in a 
strong and homogeneous static magnetic field, causing the otherwise randomly oriented 
magnetic moments of the protons, in water molecules within the body, to precess 
around the direction of the applied field. The part of the body in the homogeneous 

25 region of the magnet is then irradiated with radio-frequency (RF) energy, causing some 
of the protons to change their spin orientation. The net magnetization of the spin 
ensemble is nutated away from the direction of the applied static magnetic field by the 
applied RF energy. The component of this net magnetization orthogonal to the direction 
of the applied static magnetic field acts to induce measurable signal in a receiver coil 

30 tuned to the frequency of precession. This is the magnetic resonance (MR) signal. 

The useful RF components are those generated in at plane at 90 degrees to the 
direction of the static magnetic field. The same coil structure that generates the RF field 
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can be used to receive the MR signal or a separate receiver coil placed close to the 
patient may be used. In either case the coils are tuned to the Larmor precessional 
frequency coo where co 0 = yB 0 and y is the gyromagnetic ratio for a specific nuclide and 
B 0 is the applied static magnetic field. 
5 A desirable property of radio frequency coils for use in MR is the generation of 

homogeneous RF fields over a prescribed region. Normally this region is central to the 
coil structure for transmission resonators. A well known example of transmission 
resonators is the birdcage resonator, details of which are given by Hayes et. al. in The 
Journal of Magnetic Resonance, 63, 622 (1985) and U.S. Pat. No. 4,694,255. 

10 In some circumstances it is desirable to generate a target field over an 

asymmetric region of the coil structure, i.e., a region that is asymmetric relative to the 
mid-length point of the longitudinal axis of the coil structure. This is potentially 
advantageous for patient access, conformation of the coil structure to the local anatomy 
of the patient and for use in asymmetric magnet systems. 

1 5 One method that is known in the art for generating homogeneous fields over a 

volume that is asymmetric to the coil structure is to enclose one end of the cylindrical 
structure, a so-called 'end-cap' or dome structure (details of which are given by Meyer 
and Ballon in The Journal of Magnetic Resonance, 107, 19 (1995) and by Hayes in 
SMRM 5 th annual meeting, Montreal, Book of Abstracts, 39 (1986)). These designs 

20 were applied to structures that surrounded only the head of a patient and, by their 
nature, prevent access to the top of the head. The limited access also makes these 
structures problematic for whole-body imaging as they substantially reduce access from 
one end of the magnet. 

It is an object of this invention to provide coil structures that generate desired 

25 RF fields within certain specific, and asymmetric portions of the overall coil structure, 
preferably without substantially limiting access from one end of the structure. 
Asymmetric radio frequency coils can be used in conventional MR systems or in the 
newly developed asymmetric magnets of U.S. Patent No. 6,140,900. 

It is a particular object of the present invention to provide a general systematic 

30 method for producing a desired radio frequency field within a coil, using a full-wave, 
frequency specific technique to first define a current density on at least one cylindrical 
surface and subsequently to synthesize a coil pattern from the current density. 
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It is a further particular object of the present invention to use complex current 
densities in the full-wave, frequency specific method. 
SUMMARY OF THE INVENTION 

In one broad form, the invention in accordance with certain of its aspects 
5 provides a coil structure for a magnetic resonance device having a cylindrical space 
with open ends, the coil structure being adapted to generate a desired RF field within a 
specified portion of the cylindrical space. In accordance with the product aspects of the 
invention, this portion is asymmetrically located relative to the mid-length point of the 
longitudinal axis of the cylindrical space. 
10 In connection with another aspect, the invention provides a method for 

manufacturing a radio frequency coil structure for a MR device having a cylindrical 
space, preferably with open ends, comprising the steps of 

selecting a target region over which a transverse RF magnetic field of a 
predetermined frequency is to be applied by the coil structure, the target region being 
15 preferably asymmetrically located relative to the mid-length point of the longitudinal 
axis of the cylindrical space, 

calculating current density at the surface of the cylindrical space required to 
generate the target field at the predetermined frequency, 

synthesizing a design for the coil structure from the calculated current density in 
20 accordance with one of the methods discussed below, and 

forming a coil structure according to the synthesized design. 

Preferably, the method for calculating the current density uses a time harmonic 
method that accounts for the frequency of operation of the RF coil structure and makes 
use of a complex current density. 
25 The RF coils of the invention can be used as transmitter coils, receiver coils, or 

both transmitter and receiver coils. As discussed above, when the coil serves a 
transmitting function, it will normally be combined with a shield to reduce magnetic 
interference with external components of the magnetic resonance imaging system. To 
avoid redundancy, the following summary of the method aspects of the invention is in 
30 terms of a RF coil system which includes a main coil (corresponding to the "first 

complex current density") and a shielding coil (corresponding to the "second complex 
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current density"), it being understood that these methods can be practiced with just a 
main coil. 

In accordance with a first method aspect of the invention, which can be used 
under "mild" coil length to wavelength conditions, i.e., conditions in which the coil 
5 length is less than about one-fifth of the operating wavelength, a method for designing 
apparatus for transmitting a radio frequency field (e.g., a field having a frequency of at 
least 20 Megahertz, preferably at least 80 Megahertz), receiving a magnetic resonance 
signal, or both transmitting a radio frequency field and receiving a magnetic resonance 
signal is provided which comprises: 

1 0 (a) defining a target region (e.g., a spherical or ellipsoidal region preferably 

having a volume of at least 30 x 10 3 cm 3 and asymmetrically located relative to the 
midpoint of the RF coil, e.g., the ratio of the distance between the midpoint of the target 
region and one end of the coil to the length of the coil is less than or equal to 0.4) in 
which the radial magnetic component (e.g., B x , By, or combinations of B x and B y ) of the 

15 radio frequency field is to have desired values (e.g., substantially uniform magnitudes), 
said target region surrounding a longitudinal axis, i.e., the common longitudinal axis of 
the magnetic resonance imaging system and the RF coil; 

(b) specifying desired values for said radial magnetic component of the 
radio frequency field at a preselected set of points within the target region (e.g., a set of 

20 points distributed throughout the target region or a set of points distributed at the outer 
boundary of the target region); 

(c) defining a target surface (e.g., the shield coil surface) external to the 
apparatus on which the magnetic component of the radio frequency field is to have a 
desired value of zero at a preselected set of points on said target surface; 

25 (d) determining a first complex current density function, having real and 

imaginary parts, on a first specified cylindrical surface (i.e., the main coil surface) and a 
second complex current density, having real and imaginary parts, on a second specified 
cylindrical surface (i.e., the shield coil surface), the radius of the second specified 
cylindrical surface being greater than the radius of the first specified cylindrical surface 

30 (e.g., the radius of the second surface can be 10% greater than the radius of the first 
surface) by: 



(i) defining each of the complex current density functions as a sum of a 
series of basis functions (e.g., triangular and/or pulse or sinusoidal and/or cosinusoidal 
functions) multiplied by complex amplitude coefficients having real and imaginary 
parts; and 

(ii) determining values for the complex amplitude coefficients using an 
iterative minimization technique (e.g., a linear steepest descent or a conjugate gradient 
descent technique) applied to a first residue vector obtained by taking the difference 
between calculated field values obtained using the complex amplitude coefficients at 
the set of preselected points in the target region and the desired values at those points 
and a second residue vector equal to calculated field values obtained using the complex 
amplitude coefficients at the preselected set of points on the target surface; and 

(e) converting said first and second complex current density functions into 
sets of capacitive elements (as understood by persons skilled in the art, such capacitive 
elements will in general have some inductive and resistive properties) and sets of 
inductive elements (as understood by persons skilled in the art, such inductive elements 
will in general have some capacitive and resistive properties) located on the specified 
cylindrical surfaces by: 

(i) converting each of the first and second complex current density 
functions into a curl-free component J CU ri-free and a divergence-free component Jdiv-free 
using the relationships: 

* cur,- free = V ¥ » ^d 
^div-free = V X S , 

where \|/ and S are functions obtained from the respective first and second complex 
current density functions through the equations: 

vy=v-j 

_V 2 S = Vx J, and 

-V 2 (nS) = nVxJ, 
where n is a vector normal to the respective first and second specified cylindrical 
surfaces and J is the respective first and second complex current density functions; 
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(ii) calculating locations on the respective first and second cylindrical 
surfaces for the respective sets of capacitive elements by contouring the respective \\f 
functions; and 

(iii) calculating locations on the respective first and second cylindrical 
5 surfaces for the respective sets of inductive elements by contouring the respective 

functions n#S (v|/ and n«S are referred to herein and function as "streaming functions"). 

In accordance with a second method aspect of the invention, which can be used 
under "non-mild" coil length to wavelength conditions, i.e., conditions in which the coil 
length can be greater than about one-fifth of the operating wavelength, a method for 
1 0 designing apparatus for transmitting a radio frequency field or both transmitting a radio 
frequency field and receiving a magnetic resonance signal is provided which 
comprises: 

(a) defining a target region in which the radial magnetic component of the 
radio frequency field is to have desired values, said target region surrounding a 

1 5 longitudinal axis, i.e., the common longitudinal axis of the magnetic resonance imaging 
system and the RF coil; 

(b) specifying desired values for said radial magnetic component of the 
radio frequency field at a preselected set of points within the target region; 

(c) defining a target surface external to the apparatus on which the magnetic 
20 component of the radio frequency field is to have a desired value of zero; 

(d) determining a first complex current density function, having real and 
imaginary parts, on a first specified cylindrical surface and a second complex current 
density, having real and imaginary parts, on a second specified cylindrical surface, the 
radius of the second specified cylindrical surface being greater than the radius of the 

25 first specified cylindrical surface: 

(i) defining each of the complex current density functions as a sum of a 
series of basis functions (e.g., sinusoidal and/or cosinusoidal functions) multiplied by 
complex amplitude coefficients having real and imaginary parts; and 

(ii) determining values for the complex amplitude coefficients by 
30 simultaneously solving matrix equations of the form: 
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[Af](a c ) + [Af](o s )=S c (Eq.I) 

[Af](a c ) + [A|](o s )=B s 
CSC s 

where A 1 , Aj , A 2 , and A 2 are transformation matrices between current density 
space and magnetic field space whose components are based on time harmonic Greeks 
functions, a c and a s are vectors of the unknown complex amplitude coefficients for the 
first and second complex current density functions, respectively, B c is a vector of the 
desired values for the radial magnetic field specified in step (b), and B s is a vector 
whose values are zero, said equations being solved by: 

(1) transforming the equations into functionals that can be solved 
using a preselected regularization technique, and 

(2) solving the functionals using said regularization technique to 
obtain values for the complex amplitude coefficients; and 

(e) converting said first and second complex current density functions into 
sets of capacitive elements and sets of inductive elements located on the specified 
cylindrical surfaces (preferably, this step is performed using the methods of step (e) of 
the first method aspect of the invention; also the methods of said step (e) can be used 
independent of either the first or second method aspects of the invention to convert a 
complex current density function to sets of capacitive and inductive elements, i.e., to a 
manufacturable coil structure). 

In accordance with certain preferred embodiments of this second method aspect 
of the invention, the regularization functional is chosen so as to minimize the integral 
of the dot product of the first complex current density function with itself over the first 
specified cylindrical surface and to minimize the integral of the dot product of the 
second complex current density function with itself over the second specified 
cylindrical surface. Other regularization functionals that can be used include: (1) the 
second derivative of the complex current density functions, and (2) other functionals 
besides the dot product that are proportional to the square of the complex current 
density functions. 
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In accordance with other preferred embodiments, the complex amplitude 
coefficients are chosen so that the first and second complex current density functions 
each has zero divergence. 

The parenthetical statements set forth in connection with the summary of the 
5 first method aspect of the invention also apply to the second method aspect of the 
invention except where indicated. In connection with all of its method aspects, the 
invention also preferably includes the additional step(s) of displaying the locations of 
the sets of inductive elements on the first and second specified cylindrical surfaces 
and/or producing physical embodiments of those sets of elements. 
1 0 In accordance with certain of its product aspects, the invention provides 

apparatus for use in a magnetic resonance system for transmitting a radio frequency 
field (e.g., a field having a frequency greater than, for example, 20 megahertz, and 
preferably greater than 80 megahertz), receiving a magnetic resonance signal, or 
transmitting a radio frequency field and receiving a magnetic resonance signal, said 
1 5 apparatus and said magnetic resonance imaging system having a common longitudinal 
axis, said apparatus comprising: 

(a) a support member (e.g., a tube composed of fiberglass, TEFLON, or 
other materials, which preferably can be used as a substrate for a photolithography 
procedure for printing conductive patterns and which does not substantially absorb RF 

20 energy) which defines a bore (preferably a cylindrical bore) having first and second 

open ends which are spaced from one another along the longitudinal axis by a distance 
L (the open ends are preferably both sized to receive a patient's body part which is to be 
imaged, e.g., the head, the upper torso, the lower torso, a limb, etc.); and 

(b) a plurality of inductive elements (e.g., copper or other metallic tracks or 
25 tubes) and a plurality of capacitive elements (e.g., distributed or lumped elements) 

associated with the support member (e.g., mounted on and/or mounted in and/or 
mounted to the support member); 

wherein if used as a transmitter, the apparatus has the following characteristics 
(the "if used as a transmitter" terminology is used in defining both transmitting and 
30 receiving RF coils since by reciprocity, coils that transmit uniform radio frequency 
fields, receive radio frequency fields with a uniform weighting function): 
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(i) the apparatus produces a radio frequency field which has a radial 
magnetic component which has a spatially-varying peak magnitude whose average 
value is A r . avg ; 

(ii) the apparatus has a homogenous volume within the bore over which the 
5 spatially-varying peak magnitude of the radio frequency field has a maximum deviation 

from A r . avg which is less than or equal to 15% (preferably less than or equal to 10%); 

(iii) the homogeneous volume defines a midpoint M which is on the 
longitudinal axis, is closer to the first end than to the second end, and is spaced from 
the first end by a distance D such that the ratio D/L is less than or equal to 0.4 

10 (preferably less than or equal to 0.25); and 

(iv) at least one of said inductive elements (e.g., 1, 2, 3, ... or all) comprises a 
discrete conductor (e.g., a copper or other metallic track or tube) which follows a 
sinuous path such that during use of the apparatus current flows through a first part of 
the conductor in a first direction that has both longitudinal and azimuthal components 

1 5 and through a second part of the conductor in a second direction that has both 

longitudinal and azimuthal components, said first and second directions being different. 
An example of such a sinuous path is shown in Figure 15 discussed below. 

In certain preferred embodiments, the homogenous volume is at least 30 x 10 3 
cm 3 and/or L is at least 1 .0 meter. 
20 In other preferred embodiments, the first and second open ends have transverse 

cross-sectional areas Ai and A 2 , respectively, which satisfy the relationship: 

|Ai-A 2 |/A, <0.1. 
where Ai is preferably at least 2 x 10 3 cm 2 . 

In connection with all of the aspects of the invention, the lengths of coils and of 
25 cylindrical spaces associated therewith are defined in terms of the inductive elements 
making up the coil. In particular, for a horizontally oriented coil, the length constitutes 
the distance along the longitudinal axis from the leftmost edge of the leftmost inductive 
element to the rightmost edge of the rightmost inductive element, corresponding 
definitions applying to other orientations of the coil. 
30 BRJEF DESCRIPTION OF THE DRAWINGS 

The invention will be described by way of examples with reference to the 
drawings in which: 
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FIGURE 1 illustrates a general layout of a cylindrical radio frequency coil 
system of length L and radius a. Changing the value of D changes the degree of 
symmetry of the 'target region' relative to the coil structure; 

FIGURE 2 shows a current distribution on the surface of a cylinder calculated 
5 using a quasi-static formulation, where the vertical axis is the z-axis; 

FIGURE 3 shows plots of and J, along the length of a cylinder in the (|)=0 
and (J)=7i/2 planes respectively for a quasi-static formulation; 

FIGURES 4A and 4B show plots of J ^ and J : along the length of a cylinder in 
the (j>=0 and <J)=7i/2 planes respectively for the real solutions of Table 4 (190MHz). In 
1 0 FIGURE 4A, z o ,=20mm, while in FIGURE 4B, z O5 =50mm; 

FIGURE 5 shows the current distribution on the surface of a cylinder calculated 
from the real solution of Table 4 (z os =20mm; 190MHz); 

FIGURE 6 shows the current distribution on the surface of a cylinder calculated 
from the real solution of Table 4 (z os =50mm; 190MHz); 
1 5 FIGURE 7 shows a full conductor model of the contours of the symmetric 

current density distribution shown in FIG. 3. 

FIGURE 8 shows a current density distribution for a 320mm length coil with 
radius 100mm having a homogeneous volume 50mm from the center of the coil 
towards the coil's lower end. 
20 FIGURE 9 shows contour lines for n.S for the current density of FIG. 8. 

FIGURE 10 shows contour lines for vj/ for the current density of FIG. 8. 

FIGURE 1 1 shows the contour lines for the stream function n.S calculated 
using the second method aspect of the invention. 

FIGURE 12 shows a normalized magnetic field for a representative asymmetric 
25 coil in the transverse plane determined using the second method aspect of the invention. 
The upper curve shows the distribution in the x direction and the lower curve shows the 
distribution in the y direction. 

FIGURE 13 shows a transverse magnetic field along the z-axis for a 
representative asymmetric coil determined using the second method aspect of the 
30 invention. 

FIGURE 14 is a circuit-type layout for a representative asymmetric coil. 
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FIGURE 15 is a 3D perspective of an inductive geometry for a representative 
asymmetric coil. A partial shield is included for illustrative purposes. 

FIGURE 16 is a flow chart illustrating the first method aspect of the invention. 
FIGURE 17 is a flow chart illustrating the second method aspect of the 
5 invention. 

FIGURE 18 is a flow chart illustrating a preferred method for calculating a coil 
geometry from a current density. 

The foregoing drawings, which are incorporated in and constitute part of the 
specification, illustrate the preferred embodiments of the invention, and together with 
10 the description, serve to explain the principles of the invention. It is to be understood, 
of course, that both the drawings and the description are explanatory only and are not 
restrictive of the invention. 

DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION 

As discussed above, the present invention relates to RF coils having prescribed 

1 5 properties and to methods for designing these and other types of RF coils. Figures 1 6- 
1 8 illustrate the overall numerical procedures of the invention with reference to the 
various equations presented below. 

The methods of the invention as described below are preferably practiced on a 
digital computer system configured by suitable programming to perform the various 

20 computational steps. The programming can be done in various programming languages 
known in the art. A preferred programming language is the C language which is 
particularly well-suited to performing scientific calculations. Other languages which 
can be used include FORTRAN, BASIC, PASCAL, C**, and the like. The program can 
be embodied as an article of manufacture comprising a computer usable medium, such 

25 as a magnetic disc, an optical disc, or the like, upon which the program is encoded. 

The computer system can comprise a general purpose scientific computer and 
its associated peripherals, such as the computers and peripherals currently being 
manufactured by DIGITAL EQUIPMENT CORPORATION, IBM, HEWLETT- 
PACKARD, SUN MICROSYSTEMS, SGI or the like. For example, the numerical 

30 procedures of the invention can be implemented in C-code and performed on a work 
station. The system should include means for inputting data and means for outputting 
the results of the RF coil design both in electronic and visual form. The output can also 
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be stored on a disk drive, tape drive, or the like for further analysis and/or subsequent 
display. 

QUASI-STATIC ANALYSIS 

In outline, the object of the design process is to produce structures that generate 
5 target RF transverse fields (that is, B x or B y or combinations thereof) over a defined 
volume inside a cylinder. The value of B x (for example) is specified at a number of 
chosen points inside this region and in one embodiment, inverse methods are used to 
compute the corresponding current distribution. 

The geometry of the coil is shown in Figure 1 . Surface current densities on the 
1 0 cylinder are defined as: 

J<f) = ^ (0> *)<P + J z (<f>> z)z for p = a and -Z/2 < z < L/2 

(1) 

where r = pp + zz ; L and a are the length and radius of the cylinder. 

As is known in the art, in one formulation for the case of a symmetric coil 
1 5 structure, J(r) can be written as a summation of M sinusoidal basis functions: 



Jf(0,Z) = COS^ J Jem Sin (*,„*) 



OT=1 



(2) 



J z (^z) = sin^ £— ssLcos(k m z) 

(3) 

20 where J cm are the current coefficients and k m =(2m-\)n/L. 

However, this expression is only applicable to problems where the magnetic 
field is symmetrical about z=0. Magnetic fields are expressed in terms of magnetic 
vector potentials defined as the convolution of surface current density with the 
corresponding Green's function in the spatial domain: 

L/2 In 

25 A(r) = *&L J J G(r\r ') J(r ') <fy ' dz ' for p=a (4) 

< * 7r -A/2 0 
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where the dashed elements indicate source points and the undashed the field points. The 
corresponding generalized Green's function in cartesian coordinates is: 



5 where r and r' are the position vectors of field and source points respectively; k 0 is the 
wave-number given by kcr=2n/X 0 ; and X 0 is the free-space wavelength of the 
propagating field. 

We consider first an initial quasi-static simplification in which the wavelength 
is assumed to be very long, or the frequency near zero. Equation (5) then approaches: 



10 



G(r|r') = 



1 




(6) 



In cylindrical coordinates one solution of this is: 




(7) 



15 where r (r ' ) are again defined as the position vectors of field (source) respectively. I p 
and K p are the modified Bessel functions of 77 th order. This solution is for the case 
where observation points are internal to the cylinder, i.e., p' > p . 
The vector potentials may then be expressed as: 



L/2 In 



-in 0 



20 



(8) 



IJ2 2.7Z 



A<(r) = ^ j \G(r\r')J,(r')cosW-<P')d</>'dz 



-L/2 0 



(9) 
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L/2 2x 



A z (r) = 2*1 J $G(r\r')J;(r')Wdz' 

(10) 



1 -1J2 o 



Substituting equation (7) into the above equations and carrying out the 
integration over the cylinder's surface, we obtain the expression for (8)-(10) in the form 
of Fourier integral (series) or in terms of spectral variable k (m). The relevant magnetic 
fields are calculated from the curl of the vector potentials to be: 



B p (r) = ^ol £ e JP * \k cos kz J ,{m,k)r p (kp)K' p {ka)dk 



(11) 



B, (r) = - ±2*1 £ pe Jp* J cos fe (m> k)Ip {kp)K < p {ka)dk 

HP p=-co o 

(12) 

Writing the currents as summation of (2) and (3) and using the following 
relationship: 

B x = 5 p cos0-^sin0 

(13) 

we obtain 



B x (r) = - ^— J cos kz £ J cm vp w (*) * (ka)[kp cos 2 0/ o - cos 2#I X (kp)]dk 



(14) 

where 

r k ) - sin ^ - k «,)L/2 sin(A: + k m )L/2 
(k-k m )L/2 (k + k m )L/2 



(15) 
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The result described in equations (14) and (1 5) is known in the art as a quasi- 
static solution, being only valid in applications where the free space wavelength is 
much longer than the length of the cylinder generating the RF fields (see for example, 
H. Fujita, L. S. Petropoulos, M. A. Morich, S. M. Shvartsman, and R. W. Brown 'A 
5 hybrid inverse approach applied to the design of lumped-element RF coils,' IEEE 
Trans. Biomedical Engineering, vol-46, pp. 353-361, Mar. 1999). 

A target region is then defined as the volume within which a homogenous B x is 
required. In particular, a target field value (B xt ) is defined as the fixed magnitude of the 
desired homogenous field B x . The target region is, for example, specified as the 

10 volume inside a sphere of radius R centered at z=0. As an example, a target 

specification could be that R is set to be R « 0.7a and max(5 x )-min(^) is specified at 

less than 5%. If we set B x (r I )=B x (r 2 )=B x (r 3 )= = B x (r N )= B xt where r u r 2 , r 3 . . . .r N are 

the target field points chosen inside the target region, a set of N linear equations is 
obtained from equation (14). 

1 5 In particular, the following matrix equation is obtained: 

fi - x [ 1 U=[ r L«[^«-Li 

(16) 

where [T] is a matrix containing the values of B x at the defined observation points due 
to each of the current components; \J em ] is a column vector of the current coefficients 
20 J cl} J c2 , Jc3- . Jcm and [l] is a unity column vector. 

If the number of linear equations equals the number of unknowns (M=N), then 
the numerical solution for [«/ cm ] can be easily obtained by multiplying both sides of 
(16) by [T]-\ 

As an example of the quasi-static method, a coil of dimension: #=100 mm, 
25 Z=320 mm was designed. The specified target field in spherical coordinates was: 



Bxt — B x 


ri r 2 r 3 
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The current density coefficients up to M=3 were then calculated from the quasi- 
static expressions in equations (14) to be: 



J c i 


J c 2 




43.834274291992190 


-48.744728088378910 


36.222351074218750 



5 The corresponding resultant current distribution on the surface of the cylinder is 

shown in Figure 2 and illustrates the current paths required for the coil. In the graph of 
Figure 3, the current components J ^ and J z are plotted along the z-direction for <|>=0 

and <f>=7t/2 respectively. Contour plots (not shown) demonstrate that B x is homogeneous 
within the target region. 

10 TIME HARMONIC ANALYSIS 

In accordance with the invention, instead of using the quasi-static approach, a 
time harmonic, full wave method is used to calculate the required current density. This 
is essential in cases where the operating wavelength and the dimension of the cylinder 
are of the same order, as the prior art quasi-static method described in the above 

1 5 section is no longer valid. However, the mathematical derivation of a time-harmonic 
expression for B x in the spectral domain or an equivalent of (14) is extremely difficult. 
A more direct approach is to define B x directly in terms of the vector potentials of the 
generalized Green's function of equation (5) rather than its spectral equivalent given in 
equation (7): 



where A z are A y are the vector potentials calculated from equation (4). (Note that 
numerically, the calculation of equation (17) is carried out using a Gaussian-Quadrature 
integration routine over the prescribed 2-D cylindrical surface integral.) 

In the time-harmonic case, the wavelength or frequency of operation is taken 
25 into account and the wavenumber (k 0 ) is non-zero. This leads in many instances to the 
current density and its constituent coefficients being complex (that is, consisting of 
both real and imaginary components). In accordance with the invention, the calculated 
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current densities for the generation of specific transverse RF fields are such complex 
quantities. 

To illustrate the use of comjplex current densities, current densities were 
calculated on the same cylinder as the quasi-static example discussed above at three 
different frequencies: 190 MHz, 300 MHz and 500 MHz. The same target field was 
used, i.e., the homogeneous volume was symmetrically located. Equation (16) was 
modified to include complex current densities and a complex version of the matrix [T] 
was obtained using equation (17). The solution for [J] then follows the same 
procedures for the quasi-static example except that complex quantities are now 
included. The calculated current coefficients are shown in the following table: 



Frequency (MHz) 


Jcl 


J C 2 


Jc3 


190 


38.6895+jl.9310 


-42.7225-j2.1026 


35.0854+j 1.6609 


300 


33.4992+j5.4433 


-36.6628-j5.7392 


35.1309+j4.9677 


500 


23.8764+j 13.0168 


-25.3327-j 12.4308 


39.07675+j 14.0778 

i 



Table 1. Computed current coefficients (A/=3) at different frequencies. 



From the above table, it is clear that as the frequency increases, the imaginary 
part of the current densities becomes larger. 
BASIS FUNCTIONS FOR ASYMMETRIC SYSTEMS 

In a preferred embodiment of the invention, open-ended asymmetric structures 
are designed for a time harmonic system. For the design of asymmetric coil structures, 
surface current densities (Amp/meter) on a cylinder are defined in the same manner as 
for symmetric structures: 

J(r) = (0, z)^ + J z z)z for p=a and -1/2 <z<L/2 

(18) 

The relationship between J z and , under mild coil length/wavelength ratios, may be 
approximated as: 
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(19) 



Or alternately J z and J ^ may be decoupled and treated individually for the purposes 



of optimization. Currents on the coil are descretized along z and are represented by sub- 
domain basis functions. 

For asymmetric structures a new set of basis functions must be specified, and in 
one embodiment J z and J ^ are written as a summation of triangular T z () and pulse 
n.Q functions respectively as follows: 



(20) 



a M 

J, (<fi, z) = - sin f^-J cm {n (z - z„ + , ) -n. (z - z- )} 

*r /l = l 



(21) 



where the J cm are the current coefficients and 




(22) 




(23) 



and z m =ml : -| , z* = (m±±)/ s -| and L=L/(M + \). 
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Magnetic fields are expressed in terms of magnetic vector potentials defined as 
the convolution of surface current density with the corresponding Green's function in 
the spatial domain: 

A(r) = ZZ- f f G(r\r % )J(r % )d<f> % dz % fox p=a 

An -in o 

(24) 

As previously, the corresponding generalized Green's function in cartesian 
coordinates is: 



G(r\r<)= e 



(25) 

1 0 where r and r' are once again the position vectors of field and source points 

respectively; k 0 is the wave-number given by ko=2n/X 0 ; and A, 0 and is the free-space 
wavelength of the propagating field. 

The transverse field is written as: 

dy dz 

15 (26) 

where A z are A y are the vector potentials calculated by equation (24). 

In accordance with the invention, the homogeneous field is preferably 
asymmetrical about the z=0 plane. The target region is defined as a spherical volume of 
radius R centered at z=z os . Our first method of obtaining a set of current coefficients 
20 that produces the required homogeneity is as follows. 

Let (r/, r 2 , r 3 ... ./•#) be the N field positions chosen inside the target region and 
B x (r n ) («=1.JV) be the magnetic field at these points. The field at each target point is 
expressed in terms of M number of current coefficients as follows: 

M 

B x {r n ) = Y.C mn J cm 

25 (27) 
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where C mn represents the field contribution at r„ due to the basis function at z = z m . 

The residue (/?„) is defined to be the difference between B x (r„) and the specified 
target field B x , , i.e., 

R„{J^2-J N ) = B x {r n )-B„ n=\..N (28) 
5 A function / is defined as: 

f(J x ,J 2 ...J N ) = Y j Rl{J„J 2 J N ) (29) 

/> = ! 

In order to suppress large current variations in the solution, it is sometimes 
useful to also include the change of current amplitude as a constrained condition. In this 

10 case a term h>^aw, is added to equation (28), where m, is the difference between the 

adjacent current elements given by mr Ji - Ji+i ; w is the weight used for bringing m, to 
the same order of magnitude as R n , 

Our task is to find a set of {J, , J 2 ...J M } such that /is minimum, or 

/min = min {/(^i 5 ^ 2 - -Jm)) • There are various ways in which this function may be 
15 minimized. In one embodiment, a linear steepest descent method is used for each 

variable. An advantage of using this minimization technique is that it allows either a 
real or a complex solution to be readily obtained. 

EXAMPLES OF THE FIRST METHOD ASPECT OF THE INVENTION 

Examples of the invention, designed using the time harmonic, full wave method 
20 of the first method aspect of the invention are now detailed. 

A coil of dimension: o=100 mm, Z=320 mm was designed at different 
frequencies: 0 MHz, 190 MHz and 300 MHz. A set of target points (r„) was defined in 
spherical coordinates with respect to z OS9 where in terms of the parameters of Figure 1 , 
|z OJ | + \D\ = L/2. The radial positions of these points were r w =53.3mm; their angular 
25 positions are given in the following table: 
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1 V 


9 


Q 
t> 


1 


u 


no 
U 




HO 


no 
U 


-J 


Qf)° 


n° 
\j 


4 








n° 
u 




6 




A<\° 


7 




£0° 
OU 


g 


0° 


90° 


9 


90° 


90° 


10 


45° 


120° 


11 


45° 


135° 


12 


45° 


150° 


13 


0° 


180° 



Table 2. Angular positions of target points in spherical coordinates centered at z os . 

The specified target field was B xt =23.5 yiT. The requirement for homogeneity 
5 was specified as: max(B x )-min(B x )< 5% inside a spherical region of radius R « 0.7a for 
z O5 =20 mm and R « 0.55a for z os =50 mm. Eight current coefficients (M) were used for 
each computation. Tables 3 to 5 show the calculated current coefficients for the 
different z os values at the three operating frequencies. For the non-static cases, the real 
(Re) and imaginary (Im) components of the complex currents are shown separately. 
10 These coefficients are for use with the triangular and pulse functions discussed above. 







J C 2 


J C 3 


Jc4 


20mm 


220.68 


-188.62 


-7.60 


-42.38 


50mm 


1763.41 


-959.25 


-163.87 


-52.156 




2 os 


J C S 


Jc6 




J C 8 


20mm 


-29.62 


-36.10 


-25.10 


-62.20 


50mm 


-18.58 


-35.02 


-9.71 


-69.13 



Table 3. Computed current coefficients (M=S) for quasi-static case. 
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Zos(mm) 


Re/Im 




J C 2 


J C 3 


J C 4 


20 


Re 


226.65 


-195.19 


17.93 


-42.11 


Im 


-159.73 


+71.66 


-5.96 


-4.30 


50 


Re 


815.32 


-280.51 


-68.98 


0.22 


Im 


-505.05 


167.50 


5.94 


-2.68 




Zos(mm) 


Re/Im 


J C 5 


Jc6 


J c 7 


J C 8 


20 


Re 


-24.42 


-33.19 


-20.73 


-55.24 


Im 


+0.65 


-1.77 


-0.79 


-1.81 


50 


Re 


-33.27 


-30.49 


-9.94 


-63.67 


Im 


0.21 


0.98 


0.47 


0.44 



Table 4. Computed complex current coefficients (M=8) at 190 MHz. 

5 



Zos (mm) 


Re/Im 


Jcl 


Jc2 


J C 3 


J C 4 


20 


Re 


194.73 




-14.29 


-22.91 


Im 




117.21 


-21.10 


-1.15 


50 


Re 


634.66 


-48.97 




23.80 


Im 


0.92 


3.09 


9.78 


4.87 




Zos (mm) 


Re/Im 


J C S 


Jc6 


J c 7 


Jc8 


20 


Re 


-23.01 


-29.57 


-10.18 


-51.95 


Im 


-1.32 


-4.98 


0.61 


-10.29 


50 


Re 


-36.71 


-27.10 


-8.53 


-53.72 


Im 


4.19 


3.67 


2.50 


2.32 



Table 5. Computed complex current coefficients (M=8) at 300 MHz. 

10 The transverse magnetic field (B x ) inside the cylinder was examined for each of 

the above cases. Computed results show that the resultant B x values meet the required 
homogeneity for all the above cases and they show similar characteristics for a given 

z os . For the 1 90 MHz case, the real components of the current densities J ^ and J. are 

plotted along the z-direction in the <|>=0 and <(>=7r/2 planes on the graphs of Figure 4 A 

15 and 4B which correspond to z OJ =20mm and z OJ =50mm respectively. In both cases, J ^ 

and J z along z are of an oscillatory nature. It is seen that both components experience a 
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surge in magnitude near the positive end of the coil (i.e., the end closest to the end 
where the target region is located). 

Figures 5 and 6 show current distributions on the surface of a cylinder from the 
real solutions of Table 4 for z os equal to 20mm and 50mm, respectively. 
5 Plots of the B x field created by the currents of Figure 4 A and Figure 4B over the 

volume p < 80 mm , -160 mm < z < 160 mm (not shown) revealed a homogeneous 
region in both cases that was asymmetrical about the center of the coil and satisfied the 
criteria for homogeneity. It was also found that there was a larger volume of 
homogeneity in the z M =20mm case than in the z OJ =50mm case. 
1 0 EXAMPLES OF THE SECOND METHOD ASPECT OF THE INVENTION 

As a representative illustration of the second method aspect of the invention, the 
following sets of basis functions are used for the current density components on the 
surface of a cylinder with radius p 0 (m) and length L (m): 



q=0 p=0 n=\ m=p+\ V ^ J V ^ J 



2 ) * 



(30) 



j- =±ti 1 ^«4*+? )sJ kmZ+ Q^V (3D 



where 

20 k m =^f* (32) 

and the complex coefficients a a and c„ are to be calculated. 

This form of the basis functions was selected based on the following rationale. 
If the specified field has only one vector component that does not vary with <f> , then it 

25 can be reasoned that the term in the basis functions is unnecessary as the axis can 

be set in any direction in the xy plane. However for simplicity, <|> = 0 is set to coincide 
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with the x-axis and the term will be necessary for any specified field that has a 
y-component. 

For specified fields in MRJ applications, it is expected that the current density 
J is anti-symmetrical at 180°, so coefficients of even harmonics of <p will be zero and 
5 hence, n = 1,3,5... . Both components of J have a z dependence described with 

k m z + ^ P where /?=0,1 and m is an integer series commencing at p. The index 

p=0 describes basis functions that have J. = 0 at z = ±L/2 whereas the index p=l 
describes basis functions with a — 0 at z = ±L 1 2 . 

If the specified field is symmetric with respect to the z = 0 plane, then all 
1 0 coefficients of even m terms will equal zero and hence m can be assigned 
m = 1 ,3,5,. . . . However, all terms of m are required if the specified region is 
asymmetric with respect to the z = 0 plane as in the example coil design presented 
herein. 

The next step in the method is to calculate the coefficients of the basis functions 
15 in equations (30) and (31) such that a specified homogeneous field is generated. 
Because the time harmonic Green's functions are of a complex variable, the 
coefficients of the current density J are complex as well. If the specified field is that of 
a circularly polarized field, then the real and imaginary components would be the same 

albeit spatially separated by 90°. So if J($,z) produces Ka x intheDSV, J((|>-7r/2,z) 
20 would produce Kix y where K is just a constant. The real part of the current that would 
produce a circularly polarized field is then: 



linear =real{J((l),z)}-imag{J((l)-K/2,z)} (33) 

25 and this is the same current that would produce a linearly polarized field. Hence only 
B x k x is needed to be specified in the DSV to produce J(<j>,z) and (33) can be used to 
obtain the final current density with only a real component. 

Transmitting RF coils must have a shield to prevent eddy currents in external 
conductors and to provide a suitable RF ground. A shield may be approximated by 
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specifying B = 0 on a cylinder with a radius slightly larger than the shield radius. A 
current is approximated on the shield by the same set of basis functions as those for the 
inner coil except that the shield current now has both div-free and curl-free components 
(see below). 

5 If (a c ) and (c c ) denote column vectors of the and J z coefficients of the coil 

current (see equation (30)), (a s ) and (c s ) denote column vectors of the and J. 

coefficients of the shield current (see equation (31)), and H denotes a matrix, then the 
field in the DSV due to the coil currents is found in matrix form, as: 

10 [A](a C ) + [B](c c ) = (B-) (34) 

[C] (a s ) + [D] (c s ) = (B*) (35) 

1 5 where the summation of the vectors (B^ ) and (B* t ) results in the target field 

vector (B xt ) within the DSV. Similar matrix equations result when applying 5^=0, 

B y = 0 and B z = 0 as the condition at specified points near the outside of the surface of 
the shield cylinder: 

20 M(« c ) + [L](c c ) + [M](«^) + M(c^)=(0) (36) 

The matrix equations are usually not square because the number of points where 
the magnetic field is specified does not usually equal the number of unknown 
coefficients. The rank of the matrices are usually less than the number of unknown 
25 coefficients as well, and so a regularization method must be used for a solution. 

It should be noted that equations (34), (35), and (36) can be written in more 
compact form as Eq. (I) above. 
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REGU L ARI Z ATI ON 

The regularization method chosen as a non-limiting example was to minimize 
some functional in terms of the current density and impose the extra conditions onto the 
matrix equations such as (34) to (36). Numerous functional are available and it is 
convenient to choose a functional such that the resulting current density is easier to 
implement. The functional: 



inJj.Jrfs (37) 



min 

S 0 



results in a current density similar to one of minimum power and hence is a useful 
10 choice. When the current density equations (30) and (31) are substituted into the above 
equation (37) and differentiated to find the coefficients for a minimum: 



J oa i;tl 



S 0 



15 L^_(j.jy 5 = 0 (39) 

s 0 ° C ijkl 



the resulting matrix equation is just some constant multiplied by the identity matrix. 
This is because the basis functions are orthogonal. For a current density with no 
divergence, the wire implementation depends on the stream-function n»S, where 



q=0 p=0 n=\ m=p+] 

where n is normal to the cylindrical surface 



(40) 



, Po(Po^m a mnpq + nC ninp 

"mnpq ~ 2 ,^2,2 ^ 4 1 ) 



25 



Hence minimizing the variation in S over the surface would minimize the variation in 
its contours that ultimately determine the wire positions: 
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min j(V»V(n»S)) 2 tffe (41 A) 

So 



where V • V(n »S) gives a quantitative measure of the stream-function's variation. 
5 In this case, because the basis functions on the coil form a set with zero 

divergence, the d nmpq coefficients of S reduce to: 



d nmpq = (42) 



10 Applying the condition (42) results in a matrix with diagonal terms: 



[R] = nL 



2 ' ' v m 

J 



[/] (43) 



[R] being the regularization matrix, and [/] being the identity matrix. 
1 5 MATRIX SOLUTION 

An iterative method was used as a non-limiting example to calculate the coil 
coefficients and shield coefficients separately so that the error in each could be adjusted 
more exactly by adjusting two scalar penalty values A., and X 2 . 

Firstly, the initial coil currents were calculated using the expression: 

20 

iflf ) = {[A'f [A'] + X, [*,])"' [A'f (B x , ) (44) 

where [A'] is the combination of [A] and [B] of (34) using the relationship of equation 
(42) and [/?,] is the regularization matrix calculated by applying (37) to the coil 
25 currents. 

The field very close to the shield due solely to the coil currents is: 

(b s/c ) = [K'](af) (45) 
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where [K'] is the combination of [K] and [L] of (36). Shield currents are calculated to 
negate this field: 

5 (5) = -([M'f[M') + k 2 [R 2 ])- 1 [M'f(b s/C ) (46) 

where (sf =((a s f (c 5 ) r ), [M'] = [[M] [N]] of (36) and [/? 2 ] is the regularization 

matrix calculated by applying (37) to the shield currents. The field in the DSV due to 
the shield currents is: 

10 

(b DS v,s) = [C'](s) (47) 



where [C] = [[C] [£>]] of (35). The field in the DSV due to the coil currents is: 

15 (b DSV/c ) = [A'](af) (48) 

The error in the DSV is then: 

(b e ) - (iU-tfw)-tfw) (49) 
20 The new coil coefficient estimate is: 

{a c M ) = (a? ) + ([A'] T [A') + X, ])"' [A') 7 (b, ) (50) 

The calculation then loops back to (45) and is terminated when the difference in the 
25 error between two consecutive iterations falls below a certain predefined limit (e.g., 
0.01). The error s DSV here is defined as: 



(51) 



-29- 



and the shield error e s is defined as the error in the field produced outside the shield: 

e s =\b s ,c+[M'](s)\ (52) 

5 The penalty scalars X ] and X 2 affect the current density solution as well as determining 
the condition of the matrices to be inverted in equations (44), (46) and (50). Typically, 
the matrix problem becomes unstable (that is, where the matrix is tending towards 
being singular) below a certain value of X. 

EXAMPLES ILLUSTRATING THE ASPECTS OF THE INVENTION IN WHICH 
10 CURRENT DENSITIES ARE CONVERTED TO DISCRETE COIL STRUCTURES 
By whatever method the complex current density functions are obtained, the 
next step in the process of the overall design method is to synthesize the coil structure 
form the current density. 

In a preferred embodiment, a stream function method is employed to calculate 
1 5 the positions and dimensions of conductors such that the stream function of the current 
distribution in the conductors closely approximates the stream function of the 
calculated current density. 
QUASI-STATIC CASE 

It can be seen from equations (2) and (3) that the divergence of the quasi-static 
20 current density is zero: 

V-J = 0 (53) 
In true electromagnetic systems, the divergence is not zero but: 

VJ = -jo)p (54) 
where p is the charge density. However, for this analysis, the charge distribution is 
25 assumed to be negligible, which is true for static and low frequency systems. 

Because of zero divergence, a function S can be calculated that represents the 
magnetic flux produced from the current density as: 

VxS = J (55) 

where from equations (2) and (3): 

M c 

30 S = -ncos^-^cos^z (56) 
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where n is the radial unit vector and is normal to the cylindrical surface. Note that if J 
is a surface current density, S will only have a component perpendicular to the surface 
containing J, and hence only the magnitude of S is needed. The magnitude of S is a 
function in two dimensions. 
5 The current distribution in the conductors should be such that it produces the 

same function S as that produced by the calculated current density J. This is done by 
finding the contours of n»S and generating appropriate currents in each of the wires 
placed in the positions of the contours. If the contour 'heights' are chosen at equal 
increments between minimum and maximum, the currents in each of the wires 
10 following the contours should be equal. If the contour 'heights' are not equally spaced, 
the currents in the wires have to generated accordingly in proportion to the contour 
spacing. 

By way of an example, consider equation (56) with six real coefficients (M=6): 
c„,=47.3518, 212.791, 314.615, 687.005, -3946.81 and 5365.36. 

1 5 The values of the stream function dictates the shape and position of the contours 

and hence, is decided iteratively by investigating the electromagnetic simulations each 
contour configuration produces. For example, the contours may be taken at equal 
fractions of the stream-function's maximum value, resulting in equal current in all the 
wires. Alternatively, if the maximum is normalized to one, contours could be taken at 

20 ±0.697, ±0.394 and ±0.091 . Because the interval (0.091) between zero and n#S=0.091 
does not equal the interval (0.303) between n*S=0.091 and n»S=0.394, the current in 
the lowest 0.091 contour will have to be 0.3 (i.e., 0.091/0.303) of the current flowing in 
the 0.394 contour and 0.697 contour (these two contours should have equal current). 
The next task is to derive a workable model from these contours. 

25 Design decisions in this process include: 

1 . A magnetic symmetry plane at (j> =0 is used for the stream function in the range 
-L/2<z<L/2 and -7i<()><7t. 

2. Each contour is implemented as a conductor. Two cases were examined: one where 
the conductors are implemented as circular wires (tubes) of specified diameter and 

30 the other where the conductors are implemented as metal strips of specified width. 

The diameter or width of a conductor determines its self-inductance and hence its 
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complex impedance. It is possible to adjust the current levels in the conductors by 
adjusting the diameter or width of the conductors accordingly. 

3. The contours are connected depending on the sign of the stream function. If two 
contours are in close proximity and the stream function value has the same sign for 

5 both contours, then the contours are connected in parallel. If the contours have a 

stream function value of differing sign, the contours are connected in series. 

4. Capacitors are added at the connections of the contours connected in series and also 
at the mid-point of contours that are not connected to other conductors. However, 
there is a freedom to add capacitors wherever the designer considers suitable. The 

10 values of the capacitors are determined from an iterative approach using repeated 

frequency sweeps in the device simulation. 

5. The source points are at the extreme ends of the contour, that is, either at both ends 
or at one end: either at z=L/2 and z=-L/2 or only at z=L/2. 

A full implementation of this method results in the conductor configuration 
1 5 shown in Figure 7 for the current density defined by the six real coefficients set forth 
above, i.e., c„=47.3518, 212.791, 314.615, 687.005, -3946.81 and 5365.36. 
In a preferred embodiment, asymmetric structures are synthesized. 
When triangular basis functions are employed to approximate the current 
density, only the current-density values are provided at points within the domain of 
20 each triangular function. Thus the problem is to find the appropriate contours of a two- 
dimensional grid-like data set that will best approximate the current density that the 
data set is meant to represent. In the quasi-static problem, the current density has zero 
divergence because of the nature of the basis functions and thus the electric field at the 
surface is parallel with the current density. In the non-quasi-static case however, the 
25 electric field can have a component perpendicular to the current density and so an 
approach has to include measures for the approximation of the charge density. In 
accordance with the preferred embodiments of the invention, the curl-free and 
divergence-free components of the current density are found separately and each is 
approximated in turn. 
30 CURRENT-DENSITY COMPONENTS 

The current density J can be split into its curl-free and divergence-free 
components by writing: 
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J = VxS + V^ (57) 
From Laplace's law, it follows that 

vy=v-j (58) 

5 and the curl-free component of J is: 

J«w-/™= V V (59) 
The divergence-free component can be found by applying curl to both sides of 
equation (57) to give: 

VxVxS=VxJ 
1 0 If the divergence of S is made zero: 

VS = 0 

then: 

-V 2 S = VxJ (60) 
The divergence-free component is: 

15 J<*v-/™ =VxS (61) 

Since Jdiv^/ree is only on the surface, and to be free of singularities in the 
differentiation normal to a surface, it is best if 

nxS = 0 

Hence, (60) also becomes a Poisson equation over the surface of J: 
20 -V 2 (n.S) = n-VxJ (62) 

The task then is to solve the two Poisson equations (58) and (62) for v|/ and n*S 
and thereby using the contours of \}/ and the contours of n»S to approximate the entire 
current density J. 

Note that Jdiv-free is the current density that directly influences the magnetic 
25 field tangential to the surface of the conductor. J C urt-/ree is a measure of the quantity of 
rate of change of charge density and this quantity directly influences the electric field 
normal to the surface. In other words, idiv-free can be thought of as the inductive 
component of the circuit and S C uri-free as the capacitive component of the circuit. 

Once S and v|/ are calculated, the contours and configuration of the coils can be 
30 designed. The contours of n*S dictate the current configuration. Once n^S is found, its 
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contours and hence the coil pattern can be calculated in a similar manner to the quasi- 
static case. 

As examples of the above procedures, Figures 9 and 1 0 show the results of 
determining inductive and capacitive elements for complex current density functions 
5 obtained using the first method aspect of the invention and shown in Figure 8. The 
coil had a homogeneous region 50mm from the center towards the lower end; the 
length of the coil was 320mm and its radius was 100mm. The operating frequency was 
190MHz. 

A contour diagram of n#S is plotted in Figure 9 and demonstrates the paths for 
1 0 the conducting strips and/or wires. The contours of \\f will be demarcation curves of 
varying areas of capacitance. 

One mechanism to increase capacitance is to make an RF shield a contoured 
surface rather than just a cylindrical sheet, such that at regions of higher capacitance, 
the shield is closer to the coil than at regions where the capacitance is lower. Another 
1 5 mechanism is to add dielectric material of high permittivity at regions of higher 
capacitance compared to other regions, varying according to \\f. 

Figure 10 shows a contour plot of for the same current density as shown in 
Figure 8. The contours in this case are not current paths but are markers to indicate the 
regions of increasing charge density and hence capacitance. The region of highest 
20 capacitive effect is in the center of the contours. Note that J C uH-/ree is actually 
perpendicular to these contour lines. 

As further examples of the above procedures, Figures 1 1 through 15 show the 
results of determining inductive and capacitive elements for complex current density 
functions obtained using the second method aspect of the invention. 
25 The design objective was to produce an RF coil 20cm in diameter, 25cm in 

length with a DSV diameter of 10cm, offset by 2.5cm from the 2 = 0 plane. (The z = 0 
plane passes half-way along the length of the coil). The coil was designed at the 
frequency of 190MHz such that it could be tested in an available MRI machine. 

The second method of the invention was used with the constraint that J z = 0 at 
30 z = ±L/2. This constraint is a factor on how the resulting coil is excited because no 

current distribution flows out from the edge of the cylinder. After obtaining coefficients 
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of the current density, the function S was calculated. The contours of the stream 
function n»S for half the coil are shown in Figure 11. These are the preliminary 
patterns for conductor positions. The resulting current density was tested using a 
commercial method of moments package (e.g., the FEKO program distributed by EM 
5 Software & Systems, South Africa). The coil current density is approximated by 

Hertzian dipoles while the shield current density is ignored. Instead, a metallic shield 
approximated by triangles as per the method of moments is positioned where the shield 
should be. The normalized magnetic field of the coil in the transverse (x,y) plane is 
shown in Figure 12. The field varies within 10% over a distance of 13cm in the 

10 x-direction (upper curve) and 12cm in the y-direction (lower curve). The field variation 
along the z-axis is shown in figure 13 which shows a 10% variation over a distance of 
1 1cm shifted along the z-axis by 2.5cm. These simulated results generated by Hertzian 
dipoles approximating coil currents and a metal cylinder for the shield agree with the 
original field specification and target volume specification in the inverse program. 

1 5 The coil patterns were then converted to conductor patterns. For this case, 

2.2mm diameter wires were used for the conducting paths. Lumped elements were 
added to the model and adjusted such that the current distribution at the resonant 
frequency of 190MHz approximated the original calculated current distribution. The 
resulting coil is shown in Figure 15 with the circuit diagram shown in Figure 14. The 

20 values for the capacitors are: Cl=3.9pF, C2=4.7pF, C3=3.3pF, C4=2.7pF, C5=2.2pF, 
C6=10pF and C7 and C8 are variable capacitors. 

When the coil was constructed, it realized an unloaded Q of 139, measured from 
the 3dB down power points either side of the resonant frequency 

In all of the examples presented throughout this specification, a homogenous 

25 target field has been specified, but non-homogeneous target fields may be required in 
some circumstances and these are readily handled by the methodology disclosed herein. 

The invention has been described with particular reference to asymmetric radio 
frequency coils and a method for designing such coils, in which the section of interest 
(the DSV) can be placed at an arbitrary location within the coil. It will also be 

30 understood by those skilled in the art that various changes in form and detail may be 
made without departing from the spirit and scope of this invention. 



